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1 t Summary 


The use of probability theory to determine the effects of 
turbulent fluctuations on reaction ratv<)s in turbulent combustion 
systems is briefly reviewed. Results are presented for the effect 
of species fluctuations in particular. 'It is found that turbu- 
lent fluctuations of species act to reduce the reaction rates, in 
contrast with the temperatt^^e fluctuations previously determined 
to Increase Arrhenius reaction rate constants . ^ For the 
temperature fluctuations, a criterion is set forth for determi- 
ning if, in a given region of a turbulent flow field, the temper- 
ature can be expected to exhibit ramp-like fluctuations . Using 
the above-described results, along with results previously ob- 
tained^ ' ^ ® ® , a model is described for testing the effects of 
turbulent fluctuations of temperature and species, on reaction 

a 

rates in computer programs dealing with turbulent reactih^ flows. 
An alternative model which employs three-variable probablity 
density functions (temperature and two species) and is currently 
being formulated is discussed as well. 


2 . General Theory 

2.1 Probability Density Functions 

A pj^rameter x is said to be a continuous random vari- 
able if there exists a probability density function, p(x), 
which satisfies the following conditions'll]: 


P(x) ^ 0 (la) 

/■’■“p(x)dx = 1 (lb) 

*40 


The pdf may, of course, be defined on an interval other 
than (-“,+®), for example (0,1).* Since the pdf is defined on 
a specific interval, the functional value of p(x) is zero 
elsewhere. 

Equation (1) must be satisfied by a pdf of a one- 
dimensional continuous random variable. Probability density 
functions may be written for a multi-dimensional continuous 
random variable. A two-dimensional continuous random vari- 
able, for example, is comprised of two one-dimensiona?. con- 
tinuous random variables. A probability density function for 
a two-dimensional continuous random variable, denoted p(x,y), 
is* termed a bivariate or joint pdf. For such a pdf, the 
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conditions corresponding to equation (1) are: 


p(x,y) i 0 

_/'*^*^/‘*'*p(x,y)dx dy * 1 


(2a) 

(2b) 


The expected, or mean value of a one-dimensional con- 

f, 

tinuous random variable, x, is expressed as: 


* E(x) » p(x)dx 


(3) 


Equation (3) is also termed the first moment about the 
origin. In general, the k moment about the origin of a 
one -dimensional continuous random variable, x, is expressed 
as {ijs 


*'x,k “ (4) 

where k ~ 1,2,3,... 

For the case of k=l, equation (3) is equal to equation (4). 

The variance of one-dimensional continuous random vari- 
able, x, is expressed as: 

®x ® V(x) = _/^*(x - E(x) )^p(.x)dx (5) 
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A useful quanti.ty derived from the variance is the standard 
deviation which is the square root of the variance. Equa- 
tion (5) is the special case of the more general expression 

for the central moment, or moment about the mean, of a 

« 

one-dimensional continuous random variable, x [10]: 

pj. ^ " E(x))^p(x)dx (6) 

where k » 1,2,3,... 

Clearly, equation (6) indicates that the first central mo- 
ment (k-1) is zero. The second central moment (k-2) is 
termed the variance. For the case of k=2, equation (6) is 

■9 

equal to equation (5) . Higher central moments (k>2) are 
often used in probability theory to give further descrip- 
tions of a particular pdf under consideration. For example, 
the third central moment (k=3) , termed the skewness, is used 
to describe the symmetry or skewness of a pdf. The fourth 
central moment (k»4) , termed the kurtosls or flatness factor, 
is used to measure the "flatness" of a pdf. 

The concept of "moments*' of a one-dimensional continu- 
ous random variable may be extended to multi-dimensional 
continuous random variables. For example, the k^^ joint 
moment about the origin of a two-dimensional continuous 
random variable, (x,y), is expressed as {2]; 


where k s+r, the order o£ the moment 

r^XfZfSicct 

By comparison, equation (7) is seen to be an extension of 
equation (4). 

The joint moment about the mean of a two-dimensional 
continuous random variable, (x,y), is expressed as (2 J: 

^xy,k “ dy <B) 

where k *= s+r, the order of the moment 
s s 1,2,3, .. . 
r ~ 1,2,3,*«* 

By comparison, equation (8) is seen to be an extension of 
equation (6). 

For a specified joint pdf, it is possible to examine the 
distribution of any of the one-dimensional components of a 
multi-dimensional continuous random variable by considera- 
tion of its marginal pdf. For a joint pdf, p(x,y), the 
marginal pdf of x, for example, is given by: 

" .£^*P(x,y)dy (9) 
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The marginal distribution of x may be thought of as the 
distribution of x, with the simultaneous behavior of the 
other variable (s) suppressed. In other wordsi only the 
behavior of x is being examined. 

Utilizing the concept of a marginal pdf^ the mo- 
ment about the origin of any of the one-dimensional com- 

t 

ponents of a multi-dimensional continuous random variable 
may be expressed. For a two-dimensional continuous ran- 
dom variable (x,y), the k^^ moment about the origin of x, 
for example, is given by: • 


*‘x,k ^ 

where k - 1,2,3,... 


( 10 ) 


4- Vn 

Examination of equations (9) and (10) reveals that the k'' 
moment about the origin of x is expressed as the integral 
from -» to +®, with respect to x, of the product of x and 
its marginal pdf. 

Equation (3) may be extended to functions of a contin- 
uous random variable. In the case of a function f(x) of a 
one-dimensional random variable having a pdf p(x), the mean 
value of f (x) is: 
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Similarly, in the case of a function g(x,y) of a two- 
dimensional random variable, having the joint pdf p(x,y), 
the mean value of g(x,y) is: 


g(x,y) « P(x^y)dx dy , (12) 

The correlation coefficient, p , is a parameter 

xy 

defined for the two-dimensional continuous random variable 
(x,y) as: 


P^y " E{ Ix - E(x) Hy - E(v) ] ) 

/ V(x)V(yy (13) 

The correlation coefficient is a measure of the degree of 
linearity between x and y. Values of the correlation co- 
efficient near +1 or -1 reflect a high degree of linearity, 
while values of the correlation coefficient near zero 
indicate a lack of linearity. Positive values of the 
correlation coefficient indicate that as y increases, x 
increases. Negative values of the correlation coefficient 
indicate that y increases as x decreases. 

The numerator of equation (13) is defined as the 


covariance of x and y. The covariance is denoted by 
and expressed as$ 

“ 3t*y» « _/'^**_/'^*{x-x) (y-y)p(x,y)dx dy (X4) 
xy •• “ 

« 

By comparison of equation (10) and (8), the covariance is 
seen to be 'the second joint moment about the. mean. The 
significance of the covariance can be ascertained by con- 
sidering a two-dimensional random variable (x,y). The x 
and y are termed independent random variables if the value 
of X has no influence on the value of y (and likewise, the 

value of y has no influence on the value of x) . When x 

# 

and y are independent random* variables, the covariance if 
zero. Hence, the covariance may be considered as a minimum 
"criterion” of statistical dependence. This criterion can 
assure, at the very least, that when the covariance is not 
zero, the variables are not independent. However, no 
Statement can be made concerning independence, on this basis 
alone, if the covariance is zero [ 3 ] . 


« «4 
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2.2 Most-Likely Bivariate pdf for Two Bpettiies 

A model which accounts for the combined effects of tem- 
perature and species concentration fluctuations on the mean 
turbulent reaction rate is reviewed in this section. In this 
model, the temperature is assumed to be stwitistically indepe- 
dent of the species concentrations (Model II, ref. 4, p. 22]. 

Mathematically, this three-variable pdf for temperature 
and species is expressed as: 

p(t ,r^,rg>* f (t)g(r^,rg)., 0 i t 1 1 ^15) 

0 < r^.<l 

0 < r« < 1 

where f(t) ■ a pdf for temperature 

g(r.,r ) » a joint pdf for the concentrations 
^ ® of A and B species 

Equation (15) is a valid pdf since it satifies the following 
extension of equations (2a) and (2b) for a three-variable 

p(t,r^,r^) > 0 

o/^ o/^o'/^ dr^ dr^ - 


pdf t 
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(16a) 

(I6b) 


. .The moijb-XlHely bivariate pdf is utilized ai the 
joint pdf for the concentrations of any two species in 

this model. This pdf is selected on the basis of the excel-* 
lent agreement between the most-likely bivariate pdf, based 
on three moments, and an experimentally-measured pdf of con- 
centrations in turbulent non-reacting flow, as discussed in 
Reference 151. The three moments considered are: the first 

moment about the origin of each species concentration, the 
second joint moment about the mean (i.e., the covariance, k»2 
in equation (8)), and the third joint moment about the mean 
(k®3 in equation (8)). In addition, this pdf is selected on 
the basis of its potential for increased accuracy in the mod- 
eling of the joint pdf for species through the incorporation 
of higher moments. As an initial step in the utilization of 
the most-likely bivariate pdf for two species, the most-likely 
bivariate pdf based on the first moment about the origin of 
each species concentration and the covariance of the two con- 
centrations, is selected. The expression for the chosen form 
of the most-likely bivariate pdf for two species is: 

‘ ♦ fr • 

. 14 

gir^.rg) - 

where q «= a pKion.1 probability 

* dimensionless concentration of A 
rg ■* dimensionless concentration of B 


Th« continuous random variabXas in aquation <17) art traatad as 
passiva scalars. With this simplification, tha **q" tarm in aqua- 
tion (17) is a constant, sat aqual to unity for convanianco in 
this study. Tha valua of this constant **q” tarm is shown in 
Appendix C to have no affect on tha calculated valua of tha most- 
likely bivariate pdf and the values of tha mean quantities calcu- 
lated with it. 

The values of the constant coefficients, 
aquation (17) are obtained from the simultanious solution of the 
following constraint equations for known values of r^^, r^ and 

*^A b‘ 


0^^ 


g(rA,rg)drA - 1 

(18) 



*A 9<*A-'b>^'a ^'b ’/a 

(19) 



'b «<*'a''b’'’'a '^‘^B ■ *B 

(20) 

0^' 

0^^ 

‘■'a-'a><V'b>9<'a'‘^B>^'a '’5b - ■^A'b 

(21) 


A typical term in a reaction rate equation for species A 
would be 


...( 1 ) 

A 


W 


k(t) (r^C, ) (r, C* ) 

A A, max L> B,max. 


( 22 ) 


For the case of statistical independence between temperature and 
species concentration, the expression for the mean, turbulent reac 
tion rate 


•r(l) 

*r 

t 


w; ' « {-A (kj^t+k2)®exp [-T^/ (kj^t + kjj) 1 f (t) dt} 


(23) 


The corresponding value oC the "lavtlnar" reaction rate Is determined 
* » 

by Inserting the appropriate values of t, and r^ Into equation 
22. This yields: 




B 


- -A (kjt+kj) e*pI-V(k^t+kj)Irj,Cp 

The ratio of equations (23) and (24) Is 

• 1 . . ^ » B 


max 


<o/ (kit+kj^) exp[-T^/(kit+k2)lf (t)dt> 


(24) 


(kj^ t+k2)®exp[-T^/(kjt+k2) 1 r^ r^ 


(25) 


This ratio may be expressed as the product of a term which accounts 
for the effects of temperature fluctuations, , and a term which 
accounts for the effects of species concentration fluctuations, 

Z . These terms are expressed as follows t * 

^ a 


(kj^t+k2)^exp t-T^/ (kj^t+k2) 1 f (t)dt 


B 


(kj^t+k 2 ) exp [“T^/(kj^t+k 2 ) ] 


(26) 


Z 

r 






(27) 


The reaction rate "amplication ratio" has been determined 
4 6 8 9 

previously ' ' ' . The species "unmlxedness factor" was deter- 
mined using the numerical methods and considerations described in 

Appendices A, B, and C, Results are presented in Figures la thru 

4 

11 assuming a correlation coefficient of -0.9 . As may be seen, 
species unmlxedness predominates (low B^) when the mean concentrations 


This 


oS the speclos aro low und tui^bulent Cluotuations aitq high. 

is to be anticipated on physical grounds) that is^ as shown 

accompanying sketchy lower mean concentrations result in an 

ed probability that one or both species will be outside the 

sion volume within which reaction is possible (for the same 

* 

sity of fluctuations) . 



in the 
increase 
colli>- 
inten- 
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Also as anticipated# increasing fluctuations intensity r leads 
to increased unmixedness. 

The pdf selected for temperature is eithor the beta pdf 
or the temperature ramp’-induced pdf# based upon the criteria dis- 
cussed in Section 3. The complete model# which incorporates the 

< 

bivariate pdf formulation discussed in this Section# is detailed 
in Section 4. 


3. Development of a C^terion for the 
Selection of the Beta or Temperature Ranvp-induced PDF 

3 . 1 Experimental results from the literature 

The motivation for developing a criterion for selecting 
either the ramp pdf or beta pdf lies in observed* experimental re- 
sults. Temperature fluctuations that exhibit a ramp^like struc- 
ture have been observed experimentally under certain flow condi- 
tions# and more importantly# the flatness factor of these tempera 
ture fluctuations lies within a narrow range. 



Fiedlec^^ stvidioci a heated two-dimensional mixing layer; 
i.e,, the higher velocity stream was heated while the external 
stream had zero velocity and was at ambient temperature. A tur- 
bulent heated jet with a coaxial flowing stream was investigated 

11 12 

by Antonia« Prabhu, and Stephenson . * LaRue and Libby , and 

1 3 

Gibson, Chen, and Lin independently studied the turbulent flow 
behind a heated cylinder. Ramp-like temperature fluctuations 

* I 

were observed in dll four investigations. Even though these ex- 
periments were performed in non-reacting flows, it is assumed in 
this study that these results can be generalized to reacting flow 

as a result of the wide diversity of flow conditions under which 

13 

ramp-like temperature fluctuations were observed. Gibson et al. 

believe that this type of fluctuation is due to sharp thermal 

gradients, with the point of highest temperature* of the ramp 

occurring at either the downstream or upstream end of the turbu- 

lent/irrotational interface, depending upon the sign of the vor- 

« 

ticity of the main flow. 

Temperature signals obtained in these experiments are repro- 
duced in figure 2 from the original papers^® ' ^ . Note the 
similarity of the ramp-like structures even though they were 
obtained under different flow situations. The flatness factor 
distributions obtained by three of the researchers as they appear 

,, ,, , . g. , _ 10 , 11,12 

ed in the orrgxnal papers, are shown xn figure 3 

The abscissa' n of all graphs in figure 3 is a normalized 
radial distance, with each author using a different normalizing 


constant . 



The general trend of these curves is that the flatness factor 
has a value of approximately 3 (the Gaussian value) near the center- 
line^ decreases to a value near 2, then sharply Increases In the 
other region of the flowfield. 

This trend led Pi’e'lde/° to suggest that the value of 2 Is 
characteristic of the ’’sawtooth” appearance of the temperature 
signal. Antonia et al,^^ subsequently- agreed with Fiedler's, 
conclusion. This suggests that the flatness factor may be ’used 
as a criterion for the selection of either the beta or ramp pdf. 

4.2 An initial selection procedure 

It was pointed out in references 4 and 6 that the ramp pdf 

is only capable of generating flatness factor values between 1.0 

and 3.7. Of the four adjustable constants in the ramp pdf, the 

flatness factor is most sensitive to variations in a*. Since it 

is so sensitive, an assumption (to be tested subsequently) will be 

to assign 0* a value based upon physical grounds. Antonia and 
, . 7 

Atkinson indicate that a value of a* « 0.25 yields good agreement 
with experimental results. As a consequence of fixing 0*, varying 
ot ’ , and c over the ranges previously indicated^ generates flat- 

ness between 2,0 and 2.6. This suggests that the ramp pdf is 
applicable when F falls within this range. 

Before proceeding any further, one should ascertain whether 
this criterion agrees with experimental results. Antonia, Prabhu, 
and Stephenson (ref 11, p. 477) state that ramp-like behavior was 
observed in the region 0.6 ^ X] 1,2. The range of flatness factors 
suggested by this study, where the ramp pdf is applicable, is 2.0 
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to 2.6, If the norraullsed radial dlstaicice corresponding to this 

range of T is obtained from figure 35, one can see that 0.65 < 

— ^ 

i 

n £ 1.2. This agrees closely with the range within which Antonia 
et al,^^ observed ramp-like temperature fluctuations. This agree- 
ment supports the contention that the ramp pdf is applicable when 
2.0 £ F £ 2.6. (This implicitly justifies the assumption of con- 
stant a*.) 

An initial criterion would then be to determine if F, at a 

point in the flowfleld, lies between 2.0 and 2,6. If it ddes, . 

then the ramp pdf is applicable at that point. If it does not fall 

within that rangv^, the beta pdf is to be used. 

However, recall that the beta pdf generates the values 1.66 

< t* < 8.64. Since the beta pdf is continuous, every intermediate 

— 2 

value can be obtained by the proper selection of t and t’ . This 
implies that the beta pdf generates the values 2. 0-2. 6 for the 
flatness factor (as does the ramp pdf). As a result, the flatness 
factor cannot be used as the sole criterion due to this overlapping. 

The skewness also overlaps, but it is distinct when 2.0 £ F 
£ 2.6. For the ramp pdf, -1.57 £ S £ -1,25, and for the beta 
pdf, -0.82 £ S £ +0.82 (both when 2.0 £ F ^ 2.6). See table 1. 

The selection procedure is as follows: if the flatness fac- 

tor lies outside the range 2.0 to 2,6, then always select the beta 
pdf. If the flatness factor lies within this range, examine the 
skewness value. If -1,57 £ S £ -1,25, choose the .ramp pdfj other- 
wise select the beta pdf (see figure 5) . 
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Table 1: Summary of initial criter 

and ramp pdfs 
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The BkewneBB. alone cannot be used as a pritnaily criterion 
because the full ranges generated by the two pdf's are not dis- 
tinct . 

4 . 3 An alternative selection process 

The selection procedure outlined in the previous section is 
inconvenient from the computational point-of-view because of the 
difficulty in determining the third and fourth moments at every 
point in a flowfield. It would be preferable to use, if possible 
the first and second moments (mean and standard deviation) instead. 
This can be achieved by using the available experimental data as 
follows: 

1. On a Fvs, ri (axial distance) graph, find the 
range of r) which corresponds to 2.0 £ F <_ 2,6 

2. On a t vs. h graph, find the range of t which 
corresponds to the range of n determined xn 
step (1) . 

2 

3. Repeat step .(2) for t' , 

— 2 * 

The same ranges of t and t* result from three independent 
11 12 13 

papers ' ' after taking the different normalizing constants 

into account. The consensus is that if 

0.35 < t < 0.75 

and ^,,2 > 0.08 

then the temperature fluctuations exhibit a definite ramp-like 

character. Outside this range, the beta pdf is to be used. 

2 

It is to be anticipated that t* ^0.08 since the ramp pdf 
is applicable when the temperature fluctuations are large. It 


18 


^ ORIGINAL. .j, ■ * ^ ■ 

‘ ‘ ' ’ OF POCR Q'JA^:vr ' 

is also to be anticipated that 0.35 < F < 0,75. If t < G.35, then 
the mean tenipeirature is low, indicating a low reaction rate. This 
slow rate of reaction will not induce high temperature fluctuations. 
If F > 0.75, the mean temperature is high, implying that the magni- 
tude of the fluctuations will be a smaller percentage of the mean 
than it is at a lower temperature, effectively reducing the effect 
of the fluctuations. Hence a pdf which describes low Intensity 
fluctuations (e.g., the beta pdf) would be applicable. 


4. A Model for Determining Reaction Rates in Turbulent Reacting 
Flows 

4 

In this section, the completed Model II is set forth employing 
the elements discussed in Sections 2 and 3, as well as those discus- 
sed in previous status reports ' . 

Consider a general chemical kinetic mechanism, as shown 
below : 


k-1 


^1 ^^5 t \ ^ A 

^ ^ k-2 


I 


The "laminar" reaction rate expression for species A^ (i.e., the 
source term in the species transport equation for^A^^ would be 




2 Aj A, 


- . . .-k C + k c 

n Aj^ -n Am 


( 28 ) 
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In accordance with the model proposed here, this would be rewritten 
for a turbulent reacting flow 


Al,t 


^“r,l-2 ^ ^*'c''“t '^1^ ^“r,3-4 ^Aj^A^^”**’ 


(29) 


— ’2 

where - B^ (t, t ) is the temperature "amplification ratio" 

discussed in references 4, 6 , 8 , and 9, B j 4 ■* B . .(r.,r.,r' 

r»i-j r,i-j i y i 

2 

r! ) is the species "unmixedness " factor discussed in Section 2 

j 

herein, and k is the equilibrium constant based upon molar con- 
c 

centrations. 

In accordance with the criteria set forth in Section 3, the 
beta pdf is to be used to determine B in all regions of the tur- 
bulent flow field except when 


and 


0.35 < t 5 0.75 

, 2 > 0.08 (30) 

W 

When the conditions of eq, (30) are satisfied, the temperature ramp- 

induced pdf is to be used as detailed in ref. 9, p. 6 . 

* 

The species "unmixedness" factor, B . . , is determined from 

r,i-j 

— — ’ 2 2 

the values of c^t c ^ , cj^ and c^ obtained from appropriate trans- 
port equations. Since r. * F./c. , etc., estimates of the maxi- 

X X X ^ Tltd X 

mum molar concentrations are required. For a reactant species (fuel, 

oxidizer), c is taken to be the initial moler concentration. 

X f FI dX 

For a principal product species (eg. HO, Cp_, SiO ) , c. is taken 

» 4# X|AloX 

as the final equilibrium value determined from a preliminary equili- 
brium calculation. For an intermediate species (eg, OH, O, H) , 

c. must be estimated from a preliminary "laminar" kinetics 

i,max 

calculation using either the multi-dimensional computer program 
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of interest (eg, ^CHARNAL, SHIP) or a one-dimensional kinetics 

program (eg. ref. 14), Once F ^ , r ^ , r| and rj have been 

determined, B , . is obtained from curve-fits (or table look- 

r , i - 3 

up) of the "unmlxedness'* data presented in Section 2 herein. 

fhe abuve-described procedure ls> carried out for all of the 
terms in eq. (28) and the analogous equations for each species 
present. It should be noted that this model indicate!^ that the 
temperature fluctuations and species fluctuations result in 
opposing effects on the reaction rates. That is, B^ is always 
greater than unity, while is always less than unity. The 
behavior of these parameters (Figs, la thru 13 herein and, for 
example. Figs. 6 thru 24 of ref. 8 and Figs. 1 thru 10 of ref. 9) 
lead to the following conclusiom: 

1. In regions where the principal constituents are 


reactant species at relatively low temperature (high r^, 
r . , low t) , the principal effect will be through the in- 
creased Awhervius rate constant (beta pdf) due to the 
temperature fluctuations. 

In regions of intermediate temperature and species 
concentrations, the greatly amplified rate constant 
due to the temperature fluctuations (ramp pdf) com- 
petes with the species "unmixedness , " which may be- 
come substantial. 

In regions where combustion nears completion (high t 
and r for the principal product species) , reaction 
rates will approach their laminar values. 
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5, An Alternative Three-Variable Model for Determining Reaction 
Rates in Turbulent Reacting Flows 

In this section, a full three-variable model (temperature and 

two species), currently under development, is discussed. This model 

4 

was previously referred to as Model III . In this model, individual 
terms in the reaction rate equations have the form 


“ °l,max ^2, max *^0 “^0 '^0 *^<^^**1 P<t^X 3 ^,X 2 ) * dt dXj^ dXg 

(31) 


where the three-variable most-likely pdf is 


p(t^,Xj^,Xg) “ q * exp (Xq + Xj^t + X 2 *,+ 

\ 22 ) 

There are seven constants in this pdf, which can be evaluated by 
the following seven equations; 

/// f » / 

o 

{ 

iff dXidx^Ji - y, 

o 

f f ( ' h 

< 5 , ' .. ( 33 ) 

fff ^ ^ ^ 

^ 

f/f (X , . ? J ^ t - t ) ' /> ^^1) ^ i 

Iff ^ 

A 
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where ¥ are the moan values of species concentration and temp- 
erature and x"'"t”# x"''x ^ are the covariants of the second central 
moments. Using Newton’s method for systems > the seven non-linear 
equations can be transitormed into seven Hncer equations and the 

constants found by iteration. The equations can be rewritten as 

« • 

follows: I 

I 

b (a., X4 ‘ 

0 | 

<1 ( Ko, Xfc) • ^ , j, (K|,Xi,t)‘i'^|‘JKvclt - t * O 

-- Jx,dx».at- x;k[ =0 

« 

6 

•i (K,l,r‘'A) ’ SSV'‘r'^^(t-t)’f ‘>>‘.aHdt-x7F-0 

« 


(34) 

The equations to find the lambda's using Newton's method for the 
three-variable, most-likely pdf are as follows: 
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‘^k(, bv^ ■ ' ’ 

*^K - H. 
•^V, - b ' 


» ♦ • * 


•• 3 v;i 

^k. 


» I i 
I '« « 

i • » 


"^ki: -3 <^Vu ■ • “ 3u 


deK6K'mitM- «l>*ve 


etc. 


( 35 ) 


where is the partial derivative of equation a with respect 

Ao 

to Xq , etc . 

Each of the entries in the matricies are triple Integrals, which 
have to be numerically integrated. The determinant of the matricies 
can then be evaluated using Gauss elimination, or Gauss-Jordan 
reduction . 
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It is shown in Appendix C that for the one-variable and two- 
variable most-likely pdf's, the turbulent reaction rates are inde- 
pendent of q. As a result, the following two equations for the 
one- and two-variable pdf's can be written: 


/J k(t) p(t)dt w dfj k(t) • q p(t) dt (36) 

/J k(t) p(t,r)dt - /J /J k(t) *• q (p(t,r)) dt (37) 

» 

Noting that the Arrhenius reaction rate is the same on both sides 
ci equations (36) and (37), if the turbu)> 3 nt reaction rate is con- 
stant then the pdf with q"<l must be equal to the pdf with q»q. 
Writing the constraint equations fnr the one-variable most-likely 
pdf for q*l and q«q alternately, th.i following pattern results: 


• 1 


l 



-1 


e^* 

0 


t 




■^1 1 + vU'" 





■ e. ' 

= t 



o 

Ut 


dt ^ 


V 



<lt » 



(38) 
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By numerical solution it was shown for the one-variable most-likely 
pdf that 


V. 




(39) 


'fhin solution indicates that the constants in this pdf are changing 
in a manner which causes the pdf to remain a constant for all values 
of g. 

A similar pattern emerges when the constraint equations for 
the two-variable pdf are written with q=l and q«q alternately: 

q • s\ ^ ak= l ; V <>V ^ ^ V 

O 0 

‘L' -i ; 


(40) 

By numerical solution for the two-.variable case, it was shown that 

V 

^ I'M "I I ’ 1. ’'i. i ’"K ’’I. (41) 




This solution again shows that the constants are varyirg in a 
manner which keeps the value of the pdf a constant for any value 

of q- ■ 
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To decrease the computer time required to obtain the constants 


in the three-variable pdf^ the triple integrals in the augmented 
matrix were reduced to double integrals, where the variables are 
separable. This enables the use of the single variable Simpson's 
Rule program to evaluate the double integrals, since the integral 
of one of the variables is a standard function. The augmented ma- 
trix used to evaluate the constants in Newton's Method is composed 
of 49 triple integrals. Of the 49 integrals, 26 are duplicates, 
leaving 23 integrals that must be evaluated for each iteration., 
Because of the cross-product terms in the three-variable pdf 
equation (32) , one change of variables and two Jacobian transforms 
were utilized to change the 23 triple integrals to double inte- 
grals of the following form: 

't c- 

' where N is an integer between 0 and 5. 

When N=0 in equation ?)1, the integral over 0 is a scaled 
version of the bell-shaped distribution curve. When N*=l, this 
integral over a can be evaluated explicitly leaving a single in- 
tegral. When N=3, 4 or 5, the integral over a can be evaluated 
by parts. The resulting integrals are a function of the (Gaussian 
distribution. Details will appear in a future status report. 


27 



m 


I 


REFERENCES 


1. Mciyor, P. L., ”Xntroductory Probability and Statistical 

Applications,” Second Ed., Aidison-Wesley , Reading, MA, 

1972. 

2. ”CRC Standard Mathematical Tables,” 25th Ed., Beyer, W. H. , 

Ed., CRC Press, Boca Raton, FA,' 1980. 

3. Fry, T. C. , "Probability and its Engineering Uses,” Second 

Ed., Van Nostrand, Princeton, NJ, 1965. 

4. Antaki, P. J., Ronan, G.| Foy , E. and Chinitz, W., Semi- 

Annual Status Report (No, 3), NASA Grant NAGl-10, 2/1/81- 
7/31/81. 

5. Pope, S. B., "Probability Distribution of Scalars in Turbulent 

Shear Plow,” Proc. 2nd Symp. on Turbulent Shear Flows, 

London, 1979, pp. 3. 1-3. 6. 

6. Foy ^ E. , Renan, G. and chinitz, W. , Semi-Annual Status Rept. 

(no. 4), NASA Grant NAGl-18, 8/1/81-1/31/82. 

7. Antonia, R. A. and Atkinson, J, D. , "A Ramp Model for Tur- 

bulent Temperature Fluctuations,” Physics of Fluids , vol. 19, 
1976, pp. 1273-1278. 

% 

8. Chinj.tz, W. , Antaki, P. J. and Kassar, G. M., Semi-Annuel 

Status Rept. (No. 1), NASA Grant NAGl-18, 2/1/80-7/31/80. 

9. Kassar, G, M. Antaki, P. J. and Chinitz, W. , 'Semi-Annual Status 

Rept. (No. 2), NASA Grant NAGl-18, 8/1/80-1/31/81. 

10. Fiedler, H., "Transport of Heat Across a Plane Turbulent Mixing 

Layer," Adv. in Geophysics , vol. 19, 1974, pp. 93-109. 

11. Antonia, R. A., Prabhu, A. and Stephenson, S. E., "Conditionally 

Sampled Measurements in a Heated Turbulent Jet," J . FI . Moch . , 
.vol. 72, 1975, pp. 1956-1967. 

12. LaRue, J. 0. and Libby, Ps A., "Temperature Fluctuations in 

the Plane Turbulent Wake," Physics of Fluids, vol. 17, 

1974, pp. 1956-1967. 

13. Gibson, C. H, Chen, C. C. and Lin, S. C. , "Measurements of Tur- 

bulent Velocity and Temperature Fluctuations in the Wake of 
a Sphere." AIAA J .. Vol. 6, 1967, pp. 642-649. 

14. McLain, A. G. and Rao, C. S. R. , "A Hybrid Computer Program 

for Rapidly Solving Flowing or Static Chemical Kinetic 
Problems Involving Many Chemical Species," NASA TM X-3403, 

July 1976. 


28 


LIST OF SYMBOLS 


molar concentration of apecies i 
expected value of x 

flatness factor (or kurtosls) of a pdf 

* 

marginal pdf of x 

Arrhenius reaction rate constant 

a constant equal to T ** T . 

^ max min 

a constant equal to 

probability density function (pdf) of x 

a priori probability in the most-lilcely pdf 

normalized molar concentration equal to c./c. 

i 1 ^max 

skewness of a pdf 

activation temperature equal to activation energy 
divided by the gas constant 

a 

normalized temperature equal to (T-T . /(T -T . ) 

in i n m & X m i 

variance of x 

I 

reaction rate of species i (or a term in the reaction 
rate expression for species i) 

species "unmixedness" factor 

Arrhenius reaction rate constant amplification 
ratio 

1 

constants in the most-likely pdf 
correlation coefficient between x and y 
Gaussian mean 
Gaussian rms fluctuates 


normalized axial distance 
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Subscripts 

I 

max 

min 

t 


laminar 
maximum value 
minimum value 
turbulent 


Appendix A 


A1 Evaluation of the constant coefficients of the one- 
variable most-likely pdf, 

Al.l Numerical Procedure 

The one-variable most-likely pdf for temperature is given 
by; 

p(t) - q*exp[lg+Xj ttXgt^ ] (lA) 


where q is equal to unity, as discussed in Section 3.1. The 
constant coefficients of equation (ID), Xq,Xj,X 2 , are deter- 
mined from the simultaneous solution of the following con- 
straint equations for known values of t and t'^: 


0 

0 

0 


1 


/ p(t)dt = 1 

1 ~ 

/ t p(t)dt = t 

1 _ 

/ (t-t)p(t)dt = 



(2'A) 

(3a) 

(4A) 


Since equations (2a), (3A) and (4a) are non-linear, Newton’s 
Method [lA] for a system of non-linear equations is an 


3A , Gerald, E. F., Applied Numerical Analysis, Secoi^d Edi- 
tion, Addison-Wesley Publishing Co., Reading, MA, 1978. 
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appropriate technique of solution. This method transposes 
the original problem of solving a system of non-’linear equa- 
tions to that of solving a system of linear equations in terms 
of the "unknown” constant coefficients of equation (ID) . This 
is accomplished as follows; 

Equations (2A), ( 3 A), and (4A^ are re-written as: 


1 

= F = 0 / p(t)dt - 1 =* 0 (5A) 

1 • - 

= G ^ of t p(t)dt - t = 0 (6 a) 

H(Xorli /X2) = =. o^Nt-t) 2p(t)dt - t^ = 0 (7A) 


Let Xq=Xq, Xj-Xj and Xo=X 2 be the solution for the constant 
coefficients of equation (lA) , Let Xq^XqM 

be a point "near" the solution, (xq,Xi,X 2 )» It is 
possible to expand the functions given by equations (5A) , (6A) 
and (7 A) as a Taylor series about the point ( X g ^ 1 X | ^ 1 X 2 ^ ^ ) in 
terms of (xo-Xo'^)^ (xj-X^^^) and (x 2 -X 2 ^^). If (Xg ^ 1 x| ' 1 X 2 ^ ^ 
is "near enough" to (xq,Xj,X 2 ) it is possible to truncate the 
Taylor series after the first-derivative terms in order to 
obtain an approximate solution. This yields; 


t 


F = F + J F, (x.-xf^^) (av) 

. X . 1 1 

1=0 1 

2 

G = G + I G. (x.-XVM (9A) 

i=0 ^i ^ ^ 
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* «• 


H ** H + I H 

i-0 ^ ^ 


(lOA) 


where P. is the partial derivative of F with respect to X.. 

A • i 

3 . 

Similar ily for G and H. The functions F, G and H, and their 
associated partial derivatives are evaluated at (xJ^lxjM 



Equation (8A) , (9A) , and (lOA) form a syst'^m of linear 
equations. The "unknowns" .in these equations are the "im- 
provements" in each approximated variable: (xo-Xq^^)» 

(xi“Xi^^), (Xg-Xg'^). These equations may be solved by util= 
izing any appropriate technique for the solution of systems 
of linear equations. In the present work, Cramer's Rule is 
selected. The application of Cramer's Rule to equations (8A) 
(9A) and (lOA) yields the following expressions for the "im- 
provements" in the approximated constants of equation (lA) s 



det A . 



det A 


i = 0,1,2 
j = i+1 


(llA) 


where A - the coefficient matrix o£ the system 
of linear equations given b;;.* equation 


UIA) 




A. “ the matrix formed by replacing the 
^ elements of the column of A by 
-F, -G, -H? j * 1,2,3. 

det «* the determinant of a designated matrix 


Matrix A is expressed as: 



Solving equations (llA) for the x.'s yields: 

1 


X . 
1 



det A . 

1 

det A*' ‘ 


i - 0,1,2 
j i+1 


(12A) 


For the case of the one-variable most-likely pdf Xq,Xj,X 2 are 
unknown quantities. Hence, taken as an "ini- 

tial guess" for the solution. The resulting values for Xo,xj, 
X 2 aire "improved" solutions. Substitution of Xq,x^,X 2 for 
the values of Xq X j Xp yields "further-improved" solutions. 
The solution process continues in this manner until the de- 
sired degree of convergence is achieved. 

The expressions for the elements of matrix A are given 
in Table Section A 1.2. All integrations are performed 

using Simpson's Rule. 

*where det A is not equal to zero. 


34 


a 1.2 Expressions for the elements in matrix A 

The expressions for the elements in matrix A are ob- 
tained through the application of Liebniz’s Theorem for Dif- 
ferentiation of an Integral I2 a 1 . These expressions are 
given in Table Al. 


Table a 1. Expressions for the elements of matrix A 


Element 

Expression 


1 

q/ p(t)dt 


1 

o/ t p(t)dt 

F 

"Xj 

p(t)dt 


1 

q/ t p(t)dt 


q/ t2 p<t)dt 


q/ t^ p(t)dt 

^0 

q/ (t-t)2p(t)dt 


2A. U. S. Dept, of Commerce, Nat. Bureau of Stand., Handbook 
of Mathematical Functions with Formulas. Graphs, and 
Kathomatical Tables . AMS 55, Ed, by M. Abramowitz and I. 
Stegun, December 1972, p. 11. 
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Table Ax. (continued) 


Element 

• 

Expression 


0 / t(t-t) 2p(t)dt 



0 / (t-t) 2p(t)dt 



A1.3 Mathematical Limitations 

If the "initial guess" for the constant coefficients in 
equation (lA) is "sufficiently close" to the solution, New- 
ton’s Method converges to the solution. If the "initial guess" 
is not "sufficiently close," this method may diverge from the 
solution. Difficulty in achieving convergence is frequently 
encountered with very high and very low values of dimension- 
less mean temperature, t (e.g., t=0*l, 0.9). 

A2, Evaluation of the constant coefficients of the 
most-likely bivariate pdf 

The discussion in this section considers the most-likely 
bivariate for temperature and species. An analogous discus- 
sion applies to the most-likely bivariate pdf for two species. 
A2.1 Numerical Procedure 

The most-likely bivariate pdf for temperature and species 
is expressed as: 
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p(t,r) »* q*exp(Xo**'Xi J 


(13f;) 


wheire q is equal to unity i as discussed in Section *2. 7h'e 
constant coefficients of equation (13A) are determined from 
the simultaneous solution of the following constraint equa- 
tions for known values of t, x, and t'r*: 


/ p(t,r)dt dr ® 1 

(14 ’A) 

1 - 
/ t p(t,r)dt dr « t 

• 

(ISA) 

1 ‘ - 
/ r p(t,r)dt dr ~ r 

(16 A) 

; (t-t) (r-r)p(t,r)dt dr = t'r' 

(17-K) 


Since equations (14 a) through (17A) are non-linear, Newton's 
Method for a syf,tem of non-linear equations is utilized in the 
manner of Section Ai.l; 

Equations (14a) through (17a) are re-written ass 



= J = 

0-^ 0-^ 

p(t,r)dt dr-1 = 0 

(18A) 


= K - 

0 ^ 0 ^ 

1 ~ 
t p(t,f)dt dr-t = 

0(19A) 


= L = 

1 

o/ 0 -^ 

1 _ 
r p(t,r)dt dr-r = 

0(20A) 


= M *= 

0-^ 0 ^ 

1 

(t-t) (r-r) 



•p(t,r)dt dr-t'r* = 0 (21A) 
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Let Xo»Xo, ^3*^3 solution for the con- 

stant coefficients of equation (13‘’A) . Let Xo»Xo*l Xj^Xj^l 

X3*X3 *^ be a point “near" the solution (xq#Xj # X 2 «X 3 ) . 
It is possible to expand the functions given by equations 
(IB?\ to (23 a) as a Taylor series about the point (Xq^Ix|^ 1 
Xg^lxj*^) in terms of (xq-Xq*^), (xj-x}'^), (x 2 -X 2 ^^), 
(X 3 -Xp^). If (X 0 M X J ^ i xj * I X P ^ ) is "near enough" to 
(xo rXi f X 2 #X 3 ) f it is possible to truncate the Taylor series 
after the first-derivative terms in order to obtian an approx- 
imate solution. This yields: 


3 

j « j J 

i«0 

3 

K « K + X 
i«0 

3 

L e L + J 
i«0 

3 

M « M + I 
i~0 


J, (x.-x!^^) 
X^ 1 r 


K. (x.-x!*^) 
X . i 1 
r 


L, (x.-x!*^) 
X . XI 
1 


^X. ^*'i"^i 

X 


(22A.) 

(23A.) 

(24a) 

(25A) 


where J. is the partial derivative of J with respect to X., 

A • i 

1 

i ~ 0,1^2, 3. Similarly for K, L and M. Functions J, K, L, M 
and their associated partial derivatives are evaluated at 




/ X 2 


(l). (l) 


) 
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Equations (22a), (23K), (24^), and (2&a) forma syatem 
of linear equations. The "unknowns" in these equations are 
the "improvements" in each approximated variable: (Xo-Xo'^)V 

(X|-Xj^^)# (Xg-Xz'S# (x 3 -’x|'^). These equations are solved 
by utilizing Ci-amer’s Rule as in Section *1.1. The applica- 
tion of Cramer's Rule to equations (22*), (23h), (24A) and 
(25 a) yields the following expressions for the "improvements" 
in the approximated constants of equation (13A) : 




det B. 


JL 


det B*' ' 


i 

« 

3 


0 , 1 , 2 , 3 
i-fl 


(26A) 


where B « the coefficient matrix of the system of 
linear equations given by (26*) 

B. « the matrix formed by replacing the 
^ elements of the column of B by 
-J, -K, -L, -M; j s 1,2, 3,4. 
det = the determinant of a designated matrix 


Matrix B is expressed as: 



f where det B is not equal to zero. 
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Solving equations (26A) for the x^'s yields: 




+ 


det B 


i e 0,1, 2, 3 
j * i+1 


(27A) 


For the case of the most-likely bivariate pdf, Xq,Xi,X 2 ,X 3 
are unknown quantities. Hence taken as 

an ’’initial guess" for the true solution. The resulting 
values for xq,X},X 2 ,X 3 are "improved" solutions. Substitu- 
tion of Xo,Xj,X 2 ,Xg for the values of yields 

"further- improved" solutions. The solution process continues 
in this manner until the desired degree of convergence is 
achieved. 

The expressions for the elements in matrix B are given 
in Table 2A, Section A2.2.- These elements are initially ex- 
pressed in terms of double integrals. The large amounts of 
computer time required to solve for the constant coefficients 
of equation (13a) by utilizing this numerical procedure and 
the double-integral expressions render this solution imprac- 
tical. Use of the Monte Carlo integration technique, as op- 
posed to Simpson's Rule,, yields no appreciable decrease in 
the amounts of computer time required. In order to decrease 
this time requirement, the following method, developed during 
the course of this study, is utilized: 
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The double-integral expressions for the elements of ma- 
trix B are "reduced” to single integrals by noting, for exam- 
pie, that 

1 1 

“ 0-^ 0 ^ exp[Xo'^^l t+X^r+XatrJdt dr (283^) 

may be written as: 

1 1 

J * exp(Xo)o/ exp(Xit) [q/ exp(X2r+X 3tr)drldt (29*) 

Evaluation of the inner integral in equation (29 a) yields: 

1 exp(Xjt) 

= exp(Xo)o/ ■ - (exp(X2 + X3t)-l]dt (30^^) 

Similar manipulations arc performed with the remaining ele- 
ments of matrix B. Utilization of the single integral ex- 
pressions given in Table 3* decreases the amounts of computer 
time required and renders this solution practical. All inte- 
grations are performed with Simpson's Rule. 

D2.2 Expressions for the elements in matrix B 

The expressions for the elements in matrix B are obtained 
through the application of Liebniz's Theorem for Differentia- 
tion of an Integral [2A] . The double-integral forms of these 
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expressions are given in Table 2 a. The single-integral forms 
of these expressions are given in Table 3^. The following 
definitions are utilized in Table 2 a. 


1 cxp(X2+X3t)-l 

Xj * qS p(t,r)dr « expUo+Xjt) — 


Xj+X at 


1 


^2^0^^ P(t,r)dr * expiXQ+Xjt) 

(X2+X3t-l)exp(X2+X3t)+l 
‘ (X2+>3t)2 


Xj = p(t,r)at 


( 


(Xs+Xjt)^ - 2(X2+^st-i) - 


exp(Xo+X j t+X 2 +X 3 t) 

2 1 

exp(X2+X at) I 


(31a) 


(32a) 


(33A,) 


Table 2A. 


Element 



Double integral expressions for the 
elements for matrix B 


Expression 


1 1 

0 ^ p(t,r)dt dr 


1 1 

(,/ q/ t p(t,r)dt dr 


1 1 

q/ 0-^ ^ p(t,r)dt dr 


i 1 

q/ q/ tr p(t,r)dt dr 


1 1 

qI qS t p(t,r)dt dr 


1 1 

q/ 0 -^ p(t,r)dt dr 



Element 


Expression 



1 

q/ tr p(t,r)dt dr 


qI t^r p(t,r)dt dr 


l 

q/ r p(t,r)dt dr 


1 

qI tr p(t,r)dt dr 


q/ r^ p(t,r)dt dr 


1 

q/ tr^ p(t,r)dt dr 


1 - - 

q/ (t-t) (r-r)p(t,r)dt dr 

ft 

q/ t(t-t) (r-r) p(t,r)dt dr 


1 



q/ 0-^ r(t-t) (r-r)p(t,r)dt dr 
11 - _ 


M. 


,/ nf 


tr (t~t) (r-r) p(t,r)dt dr 
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Table 3A . Single integral expressions for the 
elements in matrix B 


Element 

Expression 


1 

of dt 



1 

of t dt 



1 

of ^2 



1 

of t X^ dt 



equal to J, 



1 

q/ t2 dt 


■ 

equal to J 

A3 


^ 3 

1 

0/ t2 X2 dt 



equal to J 

A2 



equal to J 



0/ X3 dt 


■ 

1 

0/ t X^ dt 

*• 

"xo 


+ t r J, 
AO 



+ t r J. 
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Tabic 3a. (continued) 


Element 


Expression 








L, - t L, 

- r 

t 


Xg 

X 3 A2 

X3 



A3 

0/ t2 Xj dt 

- t L, 
A3 

- r 

K, + t r J, 
As A3 

Mathematical 

Limitations 





If the "initial guess" for the constant coefficients in 
equation ( 13 A) is "sufficiently close" to the solution, New- 
ton's Method converges to the solution. If the "initial 
guess" is not "sufficiently close," this method may diverge 
from the solution. In this case, no convergence is achieved. 
Additional difficulties with convergence may arise when, in 
the denominator of X^ or X^# the value of (?i2 + X3t) is 
close to zero 


* 
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Appendix b 


The Uniqueness of the Solutions for the Constant Coefficients 
of the Most-Likely pdf 

The utilization of NGWton*^s Method for the solution 
of the unknown constant coefficients of the one'*variable 
most-likely pdf and the most»'likely bivariate pdf is 
described in detail in Appendix A. These constants 
are determined from the solution of the system of 

non-linear constraint equations for each pdf. The utilization 
of Newton's Method reduces the problem of solving a system of 
non-linear equations to that of solving a system of linear 
equations. The system of linear equations is solved by util*^ 
izing Cramer's Rule. Thus, the proof of uniqueness of the 
solutions obtained by utilizing Cramer's Rule is proof of the 
uniqueness of the solutions for the constant coefficients of 
the most-likely pdf. For simplicity, the following discussion 
demonstrates the uniqueness of the solutions for the constant 
coefficients of the one-variable most-likely pdf. An ana- 
logous discussion applies to the constant coefficients of the 
most-likely bivariate pdf. 

As shown in Appendix A , the solution to the following set 
of linear constraint equations 
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t 


P «* P + I P. 

i«0 ^ ^ 


G « G + f G (x.-x!'^ 
i«0 ^ ^ 


H •= H + I H (x.-x!'^ 
i«0 1 ^ ^ 


(8a) 


(9K) 


(lOA) 


is given by 


X. 

X 



det A . 

If 

det A 


i « 0,1,2 
j *= i+1 


( 11 ^) 


for the case where det A is not equal to zero. Equations 
(llA) are the expressions for the "improvements" in the ap- 
proximated values of the constant coefficients of the one-variable 
pdf. Cramer's Rule states that if "det A" is not equal to 
zero, then the system of linear equations given by equations 
(8A), (9A) , and (10A-) has the unique solution given by equa- 
tions (llA) , For a formal proof of this statement, the reader 

P, 

is referred to Reference 12)3 ] . Thus, the solutions for the 
constant coefficients of both the one-variable most-likely pdf 
and the most-likely bivariate pdf, obtained with the use of 
the numerical procedure outlined in Appendix A , are unique. 

‘ IB , Gerald, C. F., Applied Numerical Analysis , Second Edition, 
Addison-Wesley Publishing Co., Reading, MA, 1978. 
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Appendix d 


Variation of the Mean Turbulent Reaction Rate Constant and 
Mean Turbulent Reaction Rate with the Value of the 
a pKioKl probability, q 

The discussion in this appendix utilizes a combination 
of analytical and numerical calculations to demonstrate that 
the mean turbulent reaction rate constant obtained with the 
use of the one-variable most-likely pdf (treating temperature 
as a passive scalar) is independent of the value of the a 
pKloKi probability, q. An analogous discussion demonstrating 
that the mean turbulent reaction rate obtained with the use 
of the most-likely bivariate pdf is independent of the value 
of the a probability follows directly from that given 

for the one-variable most-likely pdf. Hence, the former dis- 
cussion is not presented. 

Let the expression for the one-variable most-likely pdf 
wherein ”q” is a constant, not equal to unity, be given by; 

f(t) = q*expUo+^ J 

Equation (IP) is utilized in Case I of this .discussion. Let 
the expression for the one-variable most-likely pdf wherein q 
is a constant, equal to unity, be given by; 



p(t) ■ explxJ+Ajt+xJt^l 


(2C) 


Equation (2c) Is utilized in Case IX of this discussion. The 
corresponding constraint equations for Cases I and II are 
given by the following: 

Case I 


1 

q/ f(t)dt * 1 

(3c‘) 

1 - 
o/ t f(t)dt t 

(4€) 

qJ (t-t)2f(t)dt = t'^ 

(5c) 

1 

q/ p(t)dt == 1 

(6P) 

1 

qJ t p(t)dt = t 

(7C) 

q/ (t-t)^p(t)dt = t 

(8C) 


The following comparison of Case I and Case II is made, utiliz 
ing the same values for t and t ' ^ for both cases, in order to 
demonstrate that the one-variable most-likely pdf, the mean 
turbulent reaction rate constants and, hence, the reaction 
rate constant amplification ratios determined for each case 
are the same. Thus: 

In order for equations (3c), (4c), and (5C) to be valid 


49 


simuXtaneously with equations (£c)> (7C) and (8C), then 

£(t) « p<t) (9c) 

4 

Equation (9C) is verified by the numerical results obtained 
with the computer program written for the procedure described 
in Appendix C. Moreover, equation (98) demonstrates that the 
values of the mean turbulent reaction rate constant determined 
for Cases 1 and XI are identical since both cases utilize 
identical expressions for the one-variable most-likely pdf. 

The numerical results of this same computer program indicate 
that this is so. 

ft 

The relationships between the constant coefficients of 
equations (IC) and (2G) may be ascertained from th? following: 
Equation (9c) indicates that 

q»explXo+Ait+Xpt^ J = exp [X q+X j t+X gt^ ] ( 108 ) 

• I 

If Xj=Xi and then 

I 4 

q«exp(Xo)exp(Xit)exp(X2t2) js 

I • • 

exp ( X 0 ) exp ( X 1 1 ) exp ( X 2 1 ^ ) 


( 118 ) 


Cancellation of "like terms” in equation (11c) yields: 

> 

I 

q*exp(Xo) * exp(Xo) (12c) 

• • 

Solving equation (12 ) for Xq yields: 

Xq s Xo - ln(q) • (l.^C) 

I i 

Equation (13 6-) as well as the fact that Xi*=Xj and X 2 *!?. axe 
verified by the numerical results obtained with the computer 
program written for the procedure outlined in Appendix a. 

Thus, knowledge of the value of the, constant a pKloh.1 
probi.bility for the most-likely pdf of a passive scalar (s) is 
not required since the values of both the pdf and the mean 
quantities calculated with it’s use are independent of "q”, 


where ”q” is a non-zero constant. 
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Pioure Z : Experimental temperature signals from Fiedler 

^ Antonia, and Gibson's original papers 
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Figure Flatness factor distribution from Antonia's original 
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Figure S : Schematic of selection process 
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